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SUMMARY 


I 


I 
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This  report  describes  the  analytical,  and  computational  aspects  of  a 
computer  program  for  predicting  Invlscld  flow  fields  and  aerodynamics  on 
realistic  reentry  configurations.  This  work  was  performed  by  members  of 
the  Mathematics  and  Engineering  Analysis  Branch  of  NSWC/WOL.  The  initial 
code  development  was  supported  by  the  Naval  Sea  Systems  Command  under  the 
Aeroballistic  Reentry  Technology  (ART)  program  with  some  of  the  fundamental 
analytical  and  numerical  work  supported  by  NSWC  Independent  Research  Funds. 
Most  of  the  final  code  development  and  documentation  was  supported  by  the 
Air  Force  Space  and  Missile  System  Organization  under  the  technical 
management  of  the  Aerospace  Corporation. 

The  authors  gratefully  acknowledge  the  efforts  of  Mr.  R.  Feldhuhn,  NSWC 
coordinator  for  the  ART  program,  who  was  responsible  for  initiating  the 
present  work  and  whose  continued  interest  and  support  throughout  the  investi- 
gation was  invaluable.  The  authors  are  also  Indebted  to  Mr.  M.  Lyons  and 
Dr.  E.  Ndefo  of  the  Aerospace  Corporation  for  several  stimulating  technical 
discussions  which  lead  to  Important  Improvements  in  the  final  code. 
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1.  Title 

A Program  for  Computing  Steady  Invlscld  Three-Dimensional  Supersonic 
Flow  on  Reentry  Vehicles;  Vol.  II:  User's  Manual 

2.  Identification 

The  following  user's  manual  Is  Intended  for  Instruction  on  how  to  use 
the  computer  program  D3CSS  (see  Ref.  1). 

3.  Configuration 

The  computer  program  was  written  for  use  on  the  CDC  6000  Series  or 
7000  Series.  The  program  is  coded  in  FORTRAN  IV  Source  language. 
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produce  a program  halt.). 

(iii.)  At  some  axial  location,  z = z^,  the  flow  field  is  completely 

known.  Tlic  plane  z = z is  referred  to  as  tlie  initial 

o 

plane  (see.  Fig.  1). 

(iv.)  The  body  shape  is  described  in  a fixed  cylindrical  coordinate 
system  sliomi  in  Fig.  1;  however,  there  are  exceptions  In 
the  case  of  a bent  nose  cone  (see  sec.  11  for  details). 
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3.  REQUIREMENTS 

1.  The  body  geometry  must  be  specified  In  a separate  subroutine 
(see.  Sec.  6). 

2.  The  flow  field  data  on  the  initial  plane  z = z must  be  supplied 

o 

on  input  TAPE3  (see.  Sec.  7). 

3.  A separate  input  card  deck  is  required  for  each  case  (see.  Sec.  7 
below) . 

4.  The  thermodynamic  properties  of  the  gas  must  be  specified  in 
subroutines  RGAS,  HRGAS,  ENTRY  ARGAS  (see  secs.  11.5  and  11.6  of 
Ref.  1). 

5.  For  a run  with  more  tangential  planes  than  25  and/or  more  radial  points 
than  20,  the  corresponding  dimension  statements  in  the  COMMON 

blocks  will  need  to  be  changed.  Also  in  the  MAIN,  D3CSS, 

the  DATA  statement:  DATA  (NCMAX=20) , (MCMAX=25)  will  need  to  be  changed. 

6.  The  program  requires  less  than  IlOK  octal  locations  of  memory  to 
run  for  a maximum  of  20  radial  points  and  25  tangential  planes. 

7.  The  program  has  the  special  recovery  routine  called  RECOVR  which 
is  supported  by  GDC  on  their  operating  systems  SCOPE  3.4  and 
KRONOS  2.1.  RECOVR  allows  the  program  to  continue  after  an  ERROR 
MODE.  Using  RECOVR  enables  the  printing  of  the  wall  pressures, 
forces  and  moments  even  if  tiiere  is  an  ERROR  MODE.  If  the  operating 
system  does  not  support  RECOVR,  then  the  RESTART  option  of  using 
TAPE15  will  have  to  be  used.  (see  Sec.  9)  or  the  Error  Mode  Update 
(see  Sec.  12)  could  be  used. 
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(i.  SPECIFICATION  OF  BODY  OEOMETRY 

A. I Tlu'  user  must  specify  the  body  geometry  in  the  subroutine  BODY. 

In  the  event  that  a body  is  to  be  computed  which  is  not  included  in  the 
version  of  BODY  in  tlie  delivery  package  (see.  Sec.  6.2),  a new  BODY 
subroutine  must  be  written.  A description  of  this  subroutine  containing 
a list  of  all  quantities  which  must  be  specitiea  is  given  in  Sec.  12.1  of 
ref.  1. 

6.2  The  version  of  BODY  contained  in  the  delivery  package  describes  a 
tamily  of  body  shapes  subject  to  the  following  criteria:** 

1.  The  baseline  body  consists  of  a nose  section  followed  by  a 
multiple  conic  with  at  most  8 conical  portions  (see  Fig.  2a). 

The  nose  section  is  either  a sphere  (tangent  to  the  first 
cone  portion)  or  a bent  sphere-cone  (see  Fig.  2e) . 

2.  The  rear  of  the  body  may  be  rounded  (see  Fig.  2b)  or  it  may 
have  an  elliptical  flare  (see  Fig.  2c).  It  may  not  have  both. 

3.  The  body  may  have  a wind  (<!^=  0),  side  (0°  < ({>*<  180°),  and 

lee  (if*=  180°)  cut  sequence.  By  a cut  sequence  at  (f  = ij)*  = PHIS 
we  mean  a sequence  of  cuts  formed  by  passing  as  many  as  four 
planes,  all  normal  to  <t)  = ({>* , through  the  body.*  (See  Fig.  2d) 

4.  A finite  span  flap  may  be  placed  on  the  second  part  of  a cut 

section; (Note:  a flap  is  modeled  as  a wedge.  At  the  end  of 

the  flap  the  wedge  is  continued  with  zero  angle  (see  Fig.  2d) 

5.  A cut  sequence  may  not  begin  after  rounding  begins  nor  on  the 
flared  section. 

The  various  parameters  required  to  define  a specific  geometry  are  read-ln 
from  cards  (see  Sec.  7). 

*A11  body  shapes  are  symmetric;  hence  a side  cut  at  (f*  implies  the  same  for 
360-(fi*, 

**A11  length  dimensions  are  assumed  non-dlmensionallzed  by  the  radius  of 
the  nose  (see  Fig.  2 e) . 
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7.  INPUT 

There  are  two  forms  of  input  to  the  program. 

7.1  Tape  Input:  The  initial  flow  field  data  required  to  begin  the 

calculation  is  on  input  TAPE3.  This  data  may  be  obtained  from  a previously 
saved  restart  file  (see,  Section  g)  or  (to  begin  the  calculation)  by 
running  a suitable  blunt  body  program  (e.g.,  that  of  ref  . 2).  The  specific 
read  statement  in  the  program  is: 

9 KtAD  (3)  NC.MC,  ATTACK,  YA»**ACH*GAMMA,HlNF»niiNF»PHlO,K,Z 

A ,NGAS,NThST ,PRX 

1 ,FN.FY,FA,MA,MY,MZ,FNZ,t-Y2,FAZ,MXZ,MYZ,*^ZZ 
? , (PHI  (M)  ,C  (M)  ,CZ  (H)  ,CPHI  (M)  ,vi=l  ,MC) 

3 , ( (R(N,H)  ,U(N,M)  ,V(N,M)  ,W(N,*3)  ,H(N,M)  ,r){N,M)  ,M=1,MC)  ,N=1,NC) 

The  variables  on  TAPE3  are  described  in  sec.  8 of  ref.  1.  We  draw  the  user's 
attention  to  the  following: 

1.  Pressures  and  densities  must  have  dimensions  consistent  with 
subroutine  RGAS;  i.e.,  lb£/ft^  and  Ib^-sec^/f t"^ , respectively. 

l/2 

Velocities  have  dimensions  of  (pressure/density)  . ATTACK  and 
YAW  are  in  degrees  but  PHI(M)  is  in  radians.  All  lengths  are  non- 
dimensionalized  by  R^  (a  characteristic  length) . Note  that  R^  must 
be  the  same  nondimensionalizing  length  as  used  in  subroutine  BODY. 

2.  The  data  points  r = R(N,M) , <p  = PHI(M)  on  this  tape  are  not  required 
to  coincide  with  the  computational  points  on  the  initial  plane  for 
the  calculation.  When  they  do  not  coincide  REZONE  must  be  used. 

It  is  however  required  that  the  array  R(N,M)  be  increasing  (i.e., 

R(N,M)  < R(N+1,M)  for  each  fixed  M)  and,  also,  the  array  PHI(M)  be 
Increasing  (i.e.,  PHI(M)  < PHI(M+1)). 
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7.2  Data  Cards  (Namelist) 

Tl.e  unformatted  data  input  to  the  program  from  cards  Is  contained 
in  three  data  decks  using  NAMELIST.  Tables  1,  2,  and  3 below  describe  the 
variables  in  each  NAMELIST  as  well  as  the  default  values  for  the  variables 


not  assigned  in  the  data  decks.  See  Appendix  B example  of  the  input 

decks  needed  to  run  the  program.  Note  that  NAMELIST  requires  that  column  1 
On  all  input  cards  be  blank. 

TABLE  1 - NAMELIST/ INPUT 1 (Input  in  MAIN) 

NA'1- LIST  / INPU^I  / K fl,  ZrnC  ,‘^ACTnP  tlTZP!-  INT  INT  ,I  ZONE, 

1 NCNtW,MCN:i'.i,IF'''  , KIK  1 N E i , I SW  SMO,  I S Wl- 00  , M 00 1 , 

2 ZMOOlO':  , 7MC01  or  , 1 3 j D , I Z J 0 , 7 0^  L 1 , 7 CEL 2 , K CE  L , KF  AO  , 

^ NJ’-lPKT,NJMKTS.'‘TA-?GET,Ta;G''T7,Trpf,pp 

if  , I3TAPT  ,KS1  A=>t^£ljj^^LCNT,1PRCFL,ISWDIF,NSGD,NSFD 


Variable  Name 

Description 

Default 

KA 

Maximum  number  of  steps  to  be  taken 

2000 

ZEND 

Last  Z value  in  the  calculation 

Undef. 

FACTOR 

Factor  is  the  CFL  factor  used  in 
computing  AZ  (c.f..  Sec.  3.6, 
ref.  1 ) 

.9 

DZPRINT 

Controls  printing  with  respect  to 
Z values;  e.g.,  the  flow  field 
variables  are  outputted  when  the 
calculation  has  marched  the  distance 
DZPRINT  along  the  Z-axis  since  the 
last  time  the  output  was  triggered. 

1.E06 

KOUT(I) 

ZPRINT(I) 

I ZONE 
NCNEW 


KOUT  and  ZPRINT  control  printing  with 
respect  to  the  number  of  steps  taken. 
KOUT(I)  is  the  number  of  steps  between 
print  outs  when  ZPRINT (I-l)i  Z < 
ZPRINT(I).  Both  vectors  are  of 
dimension  5. 

Determines  whether  or  not  to  rezone. 
IZ0NE=0,  no  rezone;  IZ0NE=1,  rezone. 

If  IZ0NE=1  then,  NCNEW  is  number  of 
points  desired  in  radial  direction 
after  re zone 


If  IZONE=l  then,  MCNEW  is  number  of 
planes  desired  in  tangential  direction 
after  rezone 


KOUT (I) =20 
for  all  I 
ZPRINT (I )= 
1.E06  for 
all  1 


0 


Undef . 


Undef . 
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Variable  Name 

Description 

I Ft: 

1 1'C  determines  order  t)f  X differ- 

encing (c.f..  Sec.  3.2,  ref.  1) 

TPC=1  tlien  foi'ward  differencing  is 

used  in  the  predictor  step  and 
backward  differencing  is  used  In 
the  corrector  step. 

n\:=()  the  differencing  is  reversed 

KIKINEC 

if  the  pressure  at  an  interior  point 
becomes  negative  tlu'  conservation  vector 
is  smoothed;  i' . f . , Sec.  4.1,  ref.  1. 

I SWSMO 

The  wall  entropv  is  extrapolated  for  the 
first  TSlv'SMO  planes  (;tero  means  no 
extrapolation);  sec  Sec.  4.2,  ref.  1. 

ISWIOD 

Determines  wliicli  version  of  wall  boundary 
condition  is  to  be  used.  ISU'MOD=0  mean.s 
MOD  0;  TSWIOD  = 3 means  MOD  3.  1 SWIOD 

must  be  0 or  3 (see  Secs.  3.4  and  10.7  of 
ref.  !) 

MOD  1 

Controls  the  use  of  second  order  accuracy 
at  tile  wall  (see.  Secs.  3.4  and  10.7 
of  ref.  1). 

MODI  = 1 means  second  order  accuracy 
will  be  used 

MOD  1 = 0 means  no  second  order  accuracy 

will  be  used  at  the  wall 

:^M0D10N 

K'hen  7.  first  becomes  larger  than 
2M0D10N,  Mi'Dl  is  set  to  1.  It  acts 
independently  of  the  pri'sent  value  of 
MOD } . 

Default 

0 


2 

0 

3 

1 


1 . F.06 
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Vari;ible_  Nanio 
ZMODIOF 

(PHTl.ID. 

I’HTJJn) 

(7,cm , 

7.CFL2) 

KCFF 

KFAC 

N.IMI’KT 

N.JMKTS 


NTARGF.T 


D e s c^r  ip t i I > n Defaiil  t 

U'lien  7.  first  bocomes  larger  than  1.E06 

7.Mi>D10F,  MODI  is  set  to  0.  It  acts 
independently  of  the  preset  value  of 
MODI. 


Defines  the  oppn  plane  Interval  measured  PHlldD=0 
in  degrees  in  whicli  JIW  SUBKOin'INF  is  PHI2.ID=0 

NOT  used  (see  Sees.  4.1  and  10.4  of 
ref.  1) 


Defines  the  opp?n  Z interval  in  which  to  ZCFhl=0 

reduce  the  CFl,  factor  (see  KFAC)  ZCFL2=0 

Determines  the  number  of  steps  after  an  0 

1‘xpansion  discontinuity  in  which  to 
reduce  the  CFi.  factor  (see  KFAC)  (see 
Sec.  4.1,  ref.  1)  KCFI.  independently  of 
(ZCFU  ,ZCF!,2) 

K'lien  tlie  CFL  factor  is  reduced,  KFAC  is  3 

the  factor  by  which  it  is  reduced, 
i.e.,  FACTOR, 'KFAC  becomes  tlie  new  CFL 
factor . 


M. IMPKT  is  tiie  maximum  number  of  steps  0 

after  an  expansion  discontinuity  for 

which  X derivatives  at  the  wall  are 
modified  (see  Sec.  4.1,  ref.  1) 

N. IMKTS  is  tlie  max.  number  of  steps  4 

after  a compression  discontinuity  for 

which  the  X derivatives  at  the  wall 
are  modified  (see  Sec.  4.1,  ref.  1) 

XTARCF.T  is  tlie  number  of  Z target  points,  0 

(Z  values  for  wh i ch  printout  is  desired,  see 
sec.  8).  There  may  he  ns  many  as  100 
t.irget  points. 


The  array  of  Z target  points.  I'ndef . 

In  the  event  of  an  t'rror  the  last  -1 

TF.RRPR  steps  will  be  printed.  If 
TF.RRPR  • 0,  then  no  printout  if  an 
error  occurs  (see  Sec.  12.4,  ref.  1). 

TSTART=0  mt'ans  do  not  restart  from  0 

TAPE  I 3 

1 START=  1 me.ins  restart  from  TAPE  1 3 
(see.  Sec.  IX). 

If  ISTART=1  then  KSTART  is  the  step  0 

number  to  restart  from  TAPI  I 3 
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Variable  Name 

Description 

Default 

ELIM 

The  error  tolerance  used  In  the 
Iterative  procedures  for  real  gas 
calculations 

.001 

LCNT 

Maximum  number  of  Iterations  for 
real  gas  Iterative  procedures 

20 

IPRCFL 

Number  of  steps  between  printouts 
of  CFL  Information 

1 

ISTOTF 

If  ISWDIF=1,  the  x differencing  Is 
switched  from  step  to 
step. 

If  ISWDIF=0,  the  x differencing  Is 
not  switched  from 
step  to  step 

0 

NSGD 

Number  of  (i>  values  to  be  Inputed 
subroutine  TRANGD  (see  sec.  10) 

in 

0 

NSFD 

Number  of  x values  to  be  Inputed 
subroutine  TRANFD  (see  sec.  10) 

in 

0 

Note:  1-Jhere  applicable  all  quantities  referring  to  axial  locations 

ai^  measured  from  the  nose  tip. 
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TABLE  2 - NA.MELIST/BODYRD  (Input  in  BODYR,  an  entry  point  of  BODY) 

NAMFL] ST/BOOyHD/NCDNE  * I WND f IEFl . ACONE f ZCONE . ZWNO t RRMU* ZFL APE ♦ 
1 THFTAFL»THFTAGL.NWtNS,NL» IFWt IPS. IFL»7W»ZS»ZLt THETAW» thetas 
? tPHlS, THETAL .HFW»hFS.HFL ♦ ZFW , ZF  S . ZF L . FA W , F A S , FAL 
3 t IBNtTMETABN,  ALNStXLVjDE  L 11  ,CENUF 


Variable  Name 

Descripti on 

Default 

NCONE 

Number  of  cone  sections,  NCONE ^ 8 

1 

IRND 

IRND=1,  if  rear  of  the  body  is  rounded; 
IRND=0.  if  rear  is  not  rounded  (see 
Fig.  2b) 

0 

I EEL 

IEFL=1,  if  rear  of  the  body  has  ellipti- 
cal flare; 

IEEL=0,  if  it  does  not  (see  Fig.  2c) 

0 

ACONE(I) 

Angle  in  degrees  of  cone  section  I,  see 
Fig.  2a  (I  < 8) 

Undef.  for 
all  I 

ZCONE  (I ) 

Z value  at  end  of  I-th  cone  section, 
see  Fig . 2a  (I  f 8) 

Undef.  I < 
1.E08,  I = 

NCONE 

NCONE 

ZRND 

If  IRNI)=1,  then  ZRND  is  the  Z value 
where  rounding  begins  (see.  Fig.  2b) 

1.E08 

RRND 

If  IRND=1,  then  RRND  is  the  radius  of 
the  round  (see  Fig.  2b) 

Undef . 

ZFLARE 

If  IEFL=1,  then  ZFLARE  is  Z value  where 
the  flare  begins  (see  Fig.  2c) 

1.E08 

THETAH., 

THETAGL 

If  IEFL=1,  THETAFL  is  the  angle  (in 
degrees)  of  growth  of  the  ((|)  = 0° , 180°) 
axis  and  THETACL  is  the  angle  (in  degrees) 
of  growth  of  the  (4>  = 90°,  270°)  axis 
(see.  Fig.  2c) 

both  undef. 

NW,NS,NL 

The  number  of  sections  in  the  wind,  side, 
and  lee  cut  sequences,  respectively. 

Note  that  NW,NS,  and  NSi  4 

0,0,0 

lEW 

1FW=1,  if  there  is  a flap  on  second 
section  of  wind  cut  sequence;  IFW=0, 
if  not. 

0 

lES 

IFS=1,  if  there  is  a flap  on  second 
section  of  side  cut  sequence; 

1FS=0,  If  not 

0 

I EL 

1FI.=  1,  if  there  is  a flap  on  second 
section  of  lee  cut  sequence; 

IFL=0,  if  not 

0 
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L 


Variable  Name 

ZW(I),ZS(I), 

7.L(1) 


THETAW(l), 
THETAS (I) , 
THETAI.(I) 

PHIS 


HFW,HFS,HFL 

ZFW,ZFS,ZFL 

FAW,FAS,FAL 


Description  Default 

Z locations  of  beginning  of  I-th  all  undef. 

section  for  wind,  side,  and  lee  cut 
sequences,  respectively.  (1^4) 

(see  Fig.  2d) 

Angles  in  degrees  of  I-th  cut  section  all  undef. 

for  wind,  side,  and  lee  cut  sequences, 
respectively.  (I f 4)  (see  Fig.  2d) 

* 

Angle  (f)  , in  degrees,  normal  to  side  90. 

cut  sequence  (see  Fig.  2d) 

Half-width  of  wind,  side,  and  lee  flaps,  all  undef. 
respectively  (see  Fig.  2d) 

Lengths  of  wind,  side,  and  lee  flaps,  all  undef. 

respectively  along  the  Z axis  (see 
Fig.  2d) 

The  angle,  in  degrees,  of  deflection,  0,0,0 

for  wind,  side,  and  lee  flaps  respectively 
(see  Fig.  2d) 


I BN 

THETABN 
XLV 
DELI  I 

CENUF 


IBN=0,  spherical  nose  section 
IBN=1,  bent  sphere-cone  nose 

Bend  angle  (see  Fig.  2e) 

Bent  nose  length  (see  Fig.  2e) 

Parameter  used  for  restarting 
on  bent  nose  (see  sec.  11) 

Parameter  used  to  test  for  axis 
shift  in  bent  nose  calc,  (see 
sec.  11) 


Undef . 
Undef. 
0 

.5 


Note  1;  It  is  not  necessary  to  set  undefined  variables  if  they  are  not 

used.  For  Instance,  if  NS=0,  THETAS  and  ZS  may  be  left  undefined. 

Note  2:  All  angles  are  in  degrees.  The  convention  used  for  determining 

whether  an  angle  is  positive  or  negative  is:  if  it  increases 

the  body  radius  with  increasing  Z it  is  positive.  If  it  decreases 
the  body  radius  witli  increasing  Z it  is  negative. 

Note  3:  Cut  sequences  must  begin  before  the  beginning  of  a rounding  or 

of  an  elliptic  flare.  Also,  the  body  cannot  be  b(\t]l  fl<Tted  and 
rounded . 

Note  4:  All  length  dimension  in  this  routine  are  nondimenslonallzed  by  the 

nose  radltis. 
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Variable 

ZREF 

AREF 

ZC 

ZO 

IPCID 
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TABLE  3 - NAMELIST/OUTRD  (Input  In  OUT) 
NAMELIST/OUTRD/ ZREF , AREF , ZC , ZO , IPCID 


Name  Description 

Reference  length  used  In  defining 

moment  coefficients  (see  Sec.  8 and 

Fig,  3).  Note,  ZREF  Is  non- 

dlmenslonallzed  by  R 
^ o 

Reference  area  used  In  definitions 
of  force  and  moment  coefficients; 
non-dlmenslonallzed  by  R^  (see 
Sec.  8 and  Fig.  3)  ° 


Z location  of  the  moment  center, 
non-dlmenslonallzed  by  Rq  (see. 
Sec.  8 and  Fig.  3) 

User  selected  origin  for  outputlng 
surface  pressure  and  aerod5Tiamlc 
coefficients  (see  Sec.  8 and 
Fig.  3) 

IPCID=0,  the  surface  pressure 
ratio  Is  printed 
IPCID=1,  the  surface  pressure 

coefficient  Is  printed 
(see  Sec,  8) 


Default 


ZENIH-ZO 


base  area  of 
the  uncut, 
unrounded, 
unflared  body 
at  z=ZEND 

0. 


0. 


0 
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7.3  Data  Cards  (Formatted) : When  the  user  chooses  to  use  a 

nonuniform  spatial  mesh  by  specifying  the  discrete  spatial  points 
(c.f.,  sec.  10),  the  required  data  is  read  from  cards  in  a formatted 
form. 

If  a nonuniform  spacing  in  the  direction  is  to  be  used,  NSGD  is 
the  number  of  points  in  the  4>  direction  (see  sec.  7.2).  The  values  of 
(KM)=PHI(M)  for  M=1,2,'»,NSGD  in  degrees  are  read  with  FORMAT (5F10. 0)  in 
subroutine  TRANGD.  Note  that  PHI(1)=0°  and  PHI (NSGD)=180°  or  360°  depending 
on  whether  the  symmetric  or  nonsyiranetrlc  problem  is  being  considered. 

If  a nonuniform  spacing  in  the  normalized  r coordinate,  x = (r-b)/(c-b), 
is  to  be  used,  NSFD  is  the  number  of  points  in  the  r direction  (see  sec.  7.2). 
The  values  of  x(N)  for  N=1 , 2 , • • • ,NSFD  are  read  with  FORMAT (5F10.0)  in 
subroutine  TRANFD.  Note  that  x(l)=0  (the  body)  and  x(NSFD)=l  (the  bow 
shock) . 
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8.  OUTPUT 

8.1  TAPE  OUTPUT : The  program  generates  two  unformatted  tapes. 

The  "WRITE"  statements  for  both  these  tapes  are  the  same,  variable  for 
variable,  as  the  READ(3)  statement  given  in  section  7.1.  The  output 
tapes  are: 

1. )  TAPE17  - This  tape  is  written  for  the  last  computational 

step  before  a normal  program  stop.  It  is  used  to  restart 
the  calculation. 

2. )  TAPE16  - This  tape  is  written  for  every  computational 

step  in  the  run.  It  is  read  in  subroutine  OUT  to  print  out 
surface  pressure  distributions  and  force  and  moment 
coefficients.  The  user  may  use  this  tape  to  interface 
with  other  programs  such  as  plotting  routines.  This  tape 
can  also  be  used  for  restarting  after  an  error  mode 
termination  of  the  program  (see  Section  9). 

8.2  Printed  Output:  The  printed  output  generated  by  the  program 

is  in  four  sections.  Example  printouts  are  given  in  Appendix  B. 

Section  1 - This  section  of  printout  is  the  heading  page  (3  copies 
are  printed).  The  heading  page  describes  the  particular  run  to 
be  made,  the  body  being  used,  and  the  various  options  selected 
for  the  calculation. 

Section  2 - In  this  section  of  output  the  flow  field  data  is 
printed  at  the  axial  locations  specified  by  the  input  parameters 

DZPRINT,  KOUT,  ZPRINT,  and  TARGETZ  (see  Table  1 of  Sec.  7.2)  Also  in 
this  section  comms^nts  and  data  associated  with  the  step  size 
determination  and  various  special  procedures  are  printed  out 


-18- 


NSWC/WOL/TR  77-32 

llircnigtuuit  the  calculation.  If  ;inerror  mode  termination  occurs, 
tills  section  contains  the  flow  field  data  for  the  last  lERRPR 
steps  when  RECOVR  is  used.  The  flow  field  data  printed  in  this 
section  contains  the  following  for  each  value  of  the  coordinate  if) 
at  the  fixed  axial  station  (note  that  (f  is  referred  to  as  ANGLE 
in  the  printout ) : 


Z axial  location  from  nose  tip 

B radius  of  body 

BZ.BPllI  derivatives  of  body  radius  with  respect  to  z and 

ip,  respectively 

C radius  of  the  bow  shock 


CZ.CPHT 


derivatives  of  the  shock  radius  with  respect  to 
z and  if,  respectively 


R radial  coordinate 

W,P,V  the  velocity  components  in  the  z,  r,  and  if  directions 

respectively  (ft /sec) 

2 

P pressure  (Ib^/ft  ) 

RllO  density  ( (lbj,-sec^) /f  t^) 

2 2 

S entropy,  free  stream  entropy  is  zero  (ft  /(sec  -°R)) 


M Mach  number 

GAMMA  effective  gamma  = ph/ (ph-p) , F in  Ref.  1 


Note:  Where  applicable,  all  quantities  are  iiondimensionalized  as  indicated 

in  Sec.  7.1. 

Section  3 - In  this  section  of  output,  either  the  surface  pressure  ratio, 

2 

p/Poo,  (when  IPCID=0)  or  the  surface  pressure  coefficient,  2(p-p^)/p^V^  , 
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(when  IPCTD=1)  are  printed  as  a function  of  ^ = z + z^^^  and  <p  where 
z^®^  = ZO  is  the  user  supplied  origin  (see  Fig.  3).  In  this  section, 
there  is  print  out  for  every  computational  point  on  the  body  surface. 

Section  A - The  aerodynamic  data  is  printed  in  two  parts,  the  first  of 

( 

which  contains  the  static  stability  coefficients  shown  in  Figure  3. 

Here  CN,  CA  and  CY  represent  the  normal,  axial  and  side  force  coefficients 
while  CMN,  CMM  and  CML  stand  for  the  pitch,  yaw  and  roll  moment  coefficients 
taken  about  x = 0,  y = 0,  z = z^.  The  variable,  XCPP  and  XCPY  are  the 
centers  of  pressure  in  the  pitch  and  yaw  planes  respectively,  measured 
from  the  origin  0^*^^  and  normalized  by  reference  length,  (see  Figure  3). 

The  printed  force  and  moment  coefficients  are  defined  by: 
force  coeff  = (force)/q 
moment  coeff  = (moment) /q  A^^^. 

Here  q = and  A^^^  and  are  supplied  by  the  user  in  NAMELIST  OUTRD. 

The  above  quantities  are  printed  at  ZEND  and  the  axial  locations  specified 
by  TARGETZ.  In  the  print  out,  the  axial  location  are  referred  to  the 
origin  see  Fig.  3. 

In  the  second  part  of  this  section,  the  z derivatives  of  the  force 

and  moment  coefficients  are  printed  as  a function  of  ^ = z + z^^^  for 

every  axial  step  in  the  calculation.  The  notation  in  this  second  part 

follows  that  of  the  first  part;  e.g.,  CNZ  = etc. 

0 z 

IMPORTANT  NOTE:  The  user  selected  origin  is  used  only  in  the  printout  of 

the  aerodynamic  data.  The  code  always  operates  in  the  coordinate  system 
with  origin  at  0 (see  Fig.  3).  Therefore,  the  body  geometry  and  all 
program  controls  must  be  referred  to  the  origin  0. 
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9 . RESTART  AND  REZONE  CAPABILITIES 

If  desired,  the  program  may  be  stopped  before  the  end  of  the 
body  and  restarted  and/or  rezoned.  This  option  is  used,  for  example, 
to  change  options  to  handle  a different  characteristic  of  the  flow. 

Another  possible  reason  for  restarting  is  to  change  the  mesh.  This  is 
done  through  the  rezone  option  (see  IZONE  in  the  input  data  section). 

Tlie  rezone  operation  interpolates  (linearly)  to  determine  now  flow  field 
variables;  c.f..  Sec.  11.3  of  ref.  1.  It  is  inadvisable  to  rezone  after 
a discontinuity  in  body  slope  has  been  encountered  because  the  large 
gradients  in  the  flow  field  which  result  will  make  the  interpolated  data 
erroneous . 

The  program  allows  the  user  two  options  for  restarting.  The  steps 
for  the  first  option  are  as  follows  (see  Fig.  Al: 

1.  After  the  first  run,  save  TAPE17. 

2.  For  the  next  run,  use  TAPE17  from  the  previous  run  as  the 

input  tape,  TAPE3.  Again  if  desired,  TAPE17  can  be  saved. 

3.  For  further  restart  runs,  repeat  step  2. 

The  steps  for  the  second  option  are  as  follows  (see  Fig.  5): 

1.  For  the  first  run,  REQUEST,  TAPE16  as  an  output  tape. 

2.  For  the  next  run,  use  TAPE  16  from  the  previous  run  as  the 

input  tape,  TAPE15.  Read  in  1START=1  and  KSTART  as  the 
step  number  at  which  to  restart  (see  ISTART  and  KSTART  in 
the  input  data  section).  If  desired,  TAPE16  may  again 
be  requested  as  an  output  tape. 

3.  For  further  restart  runs,  repeat  step  2. 

The  program  also  has  provisions  for  the  automatic  restarting  and 

rezonlng  required  after  an  axis  shift  In  bent  cone  calculations  (see  sec.  11). 

Note : In  bent  cone  calculations,  the  user  can  restart  and  rezone  without  an 

axis  shift  using  the  above  options,  however.  In  these  cases  core  must  be 
zeroed  before  the  restart  run. 
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' ""input  ^ 

DATA  CARDS 

TAPE 17 


PRINTED 

OUTPUT 


Fig.  4,  Restart  Using  First  Option 
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Fig.  5,  Restart  Using  Second  Option 
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10.  MESH  CLUSTERING 

The  program  provides  two  methods  for  clustering  the  computational 
points  in  the  shock  layer  in  either  the  r or  4i  or  both  directions.  The 
first  uses  analytic  mapping  functions,  the  other  allows  the  user  to 
directly  input  the  desired  mesh  points  in  the  shock  laver.  Both  methods 
are  discussed  in  section  4.3  of  Ref.  1.  Note  that  the  clustering  feature 
is  contained  entirely  in  the  subroutines  TRANE  and  TRANG  (see  sec.  10,5 
and  10.6  of  Ref.  1 for  a description  of  these  routines). 

10.1  Analytic  Mesh  Clustering  Functions:  At  the  present  time, 

we  have  not  sufflcientlv  studied  the  use  of  analytic  mesh  clustering 
functions  and  therefore  cannot  make  specific  recommendations.  For  the 
convenience  of  the  user  who  may  wish  to  experiment  with  the  use  of  this 
feature,  we  Include  here  Instructions  for  incorporating  mesh  clustering 
functions  into  the  program.  The  first  step  is  to  select  the  mapping 
function  f(X,Y,Z)  or  g(Y,Z),  or  both,  which  produce  the  desired  clustering 
and  which  are  admissible.  The  restrictions  on  the  functions  f and  g and 
an  example  are  given  in  .Sec.  4.3  of  ref.  1.  If  the  selected  function  f 
is  not  frX  then  the  delivery  version  of  TRANF  must  be  modified.  If 
the  function  g selected  is  not  giY  then  the  delivery  version  of  TRANG 
must  be  modified.  For  easy  reference,  listings  of  the  delivery  versions 
of  TRANF  and  TRANG  are  given  in  Figures  6 and  7,  respectively. 
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To  modify  TRANF , FORTRAN  expressions  must  he  written  wiiich  define 
the  following  variables  as  functions  of  X(N),  Y(M),  and  Z: 

SX  = f , SFX  = f^  , SFY  = f^  , SFZ  = f^ 

SFXX  = f^^  , SFYX  = f^_^  , SFZY  = f^^^ 

TRANF  can  tl<en  be  modified  using  the  following  UPDATE  IDENT: 

*IDENT  TRANFl 

*DEEETE  TRANF . 28 .TRANF . 29 

I me ‘it  FORTRAN  cards  ^or  computing  sfx.sfy,  etc.  which 
do  not  depend  on  x(N) 

*DELETE  TRANF. 40 

Insert  FORTRiVN  cards  ^or  computing  SX  and  SFX.SFY,  etc. 
which  depend  on  x(N) 

*DELETE  TR.\NF.61 

3010  FORMAT(llX,*Irfcuti J^tcatu'u  Statement  ){or  new  transformation 
function  f*) 

To  modify  TRANG,  FORTRAN  expressions  must  be  written  which  define 
the  following  variables  as  functions  of  YY(=Y)  and  Z 

SG  = g , SGY  = g^,  , SGZ  = g^^  . 

SGY^  = g^Y  , SGYZ  = g^j, 

TRANG  can  be  modified  using  the  following  UPDATE  lOENT: 


* IDENT  TRANG 1 
*DE1,ETE  TRANG.  38 

SG= 

SGY= 

SGZ=  > 

SGYY= 

SGYZ=  , 


(Expressed  as  functions  of  YY  and  z 


*DEI,ETE  TRANG.  54 

3010  FORMAT(llX,*I(te»it{,'((Crtt<o»[  Statement  for  new  transformation  g*) 
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SOaiJOUTINt  THANF(M«J,n 

c 

C TK4NK  OEFlMtS  OUANTlIItb  AbSOCIflltD  KITH  THt  CLUSTEHIN6 

C ruANSFOHUATION  IN  THE  H OIHECIION  (SEE  STATEMENTS 

C 1-8  RELOmI.  TmE  CLUSTEHING  THANSFOUMArinN  IS  ASSUMED 

C IN  THE  FO~m 

C SX  = SF(X.T,2)  whEME  S*=(»-B(7fPHin/(C(7.PMI)-BlZtPHin 

C The  user  can  SPFCI'T  the  function  SF(*.Y.2). 

C The  user  must  define  as  functions  of  (X.r.z)  the  FOLLOWING 

C SA.  SFX.  SFTt  SFZ.  Sfxx,  SFZX.  SFVX 

C SEE  USERS  manual  FOR  RESTmICHOnS  aNO  INSTRUCTIONS 

C J=1,Z.3  IS  A line  index  FOR  THANF  (JUANTITIES 

C 1=1.2  IS  A LINE  INDEX  FOR  IhanF  QUANTITIES 

C 

common  NC.mC.K.PINF.DINF.PHIO. IOYAw.PI.RAO 
COMMON  YZOl  .YPhI  (3)  .C  (2bi  .CZ  (251  .CRH  I ( 25 ) . R ( 20 . 25 1 
COMMON  0(2U.2S).P(20.25).0(2D525).V(20,?S).i«(20.2S!.AS(3(20.25I 
common  CU(4.20.2S) .CUP(A,20.25) 
c •••  END  OF  blank  common 

COMMON  /CTRANF/  NSFD.SFD(20! .SFxO(20) <5FXXO(201 

common  /CHOOY/  Z.BZZ.BPhPHI .BZPHI .tanco.oelz 

1 ,PhI (2SI ,H(2S) .nZ(2S) vHPhI (25) .COSPh! (25) .SINPMI (25) 
common  /BLK02/  ThETA.OY.TGa  ( 3)  , TG5  ( 3)  t IGF.  ( 3) 

J .x(20).XZ(20.2).XM(20.2).XPMl(2n.2).V!?S) 

2 ,TF»(20t?>.TFA(20.2».TF7(20-.  2) 

C 

CZMrCZ(M)  i BM  = B(M)  » t3ZMrHZ(M)  18PH I M = BPH I ( M ) 

CMb  = C(M)-HM  s t)ZMCZ  = BZM  - C ZM  % B?MCP  = BPh1M  -CPHI  (M) 
YZJ=YZ(J)  S YPHI JrYPHI ( J)  S TG5J=TG6(J) 


THIS  routine  is  version  1 of  twanf  corresponding  either 
TO  no  clustering.  I.E..  SF(X.Y.Z)=X 
OH  THE  USER  HAS  READ  IN  TmE  SF ( X ) DATA  POINTS 

SEX=1,0  % SFXX=0.0  S SFTrO.u 

SFZx  = 0.0  S SFYX  = 0.0  S S(Z  = 0,0 

DO  100  N=I ,NC 
IF  (NSFD  .EO.  0)  go  to  25 

The  user  read  in  the  sf(xi  data  points 

SX=SFD(N)  t SFX=SFX0(N1  i SFXXsSFXVO(N) 

GO  TO  50 

CORRESPON05  TO  NO  clustering 
25  SX=X(N) 

C •••  The  following  STATEMENTS  SHOULD  APREAR  IN  ALL  VERSIONS  ••• 
50  5X1=SX-1. 

FX=l ,/SFX 

F (H  T = -SFY»Fa«YPHI J 
FZ  = -FX« (SF  Z»SF  Y»TZJ) 

1 H (N.M)  =r,M.bX»CMR 

2 XM (N. I ) =TXR^FX/CMb 

3 XZ IN. 1 ) -FZ -TAR* (SXl •BZM-SX»CZM1 

A XPHI(N.I)rFrHD»TXR« ( 5 X 1 »HRH I M - b x • CPH 1 (Ml ) 

5 TF  4 ( N . I ) =SF  XX/SF  X 

I TF h (N. I ) =TG5J  » (5F ZX»SF YX*YZJ) /SFX- PZMCZ/CmB 
H TF / (N. I ) =SF  YX'YPHI J/SFX-BPMCP/CMb 
100  CONTINUE 
RETURN 

entry  TRANFii 
IVERSON=l 

WRITE  (F..3000)  IVERSON 
IF  (NSFO  .EO.  01  WRITE  (6.30101 
IF  (NSFO  .NE.  01  WRITE  (5.3020) 

3000  format ( 1H0.20X. •PROGRAM  TRANF • . 6X . • VERS 1 ON» . I ♦ 1 
3010  format ( I IX. •EQUAL  SPACING  IN  RA0IA(  OIRFCTION^) 

3020  format ( 1 IX .«SF ( X ) WAS  REAO  IN  AS  DATA  POINTS*) 

RETURN 

END 
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subroutine  THANGITYtMtJ) 

Trang 

2 

c 

trang 

3 

c 

IHANG  defines  OUANTITIES  associated  »ITH  the  clustehing 

trang 

4 

c 

ThaNSFORmat ION  IN  THE  PHI  direction  (SEE  STATEpENTS 

trang 

5 

c 

1-G  HELOti).  THE  CLUSTERING  TRANSFORMATION  IS  ASSUMED 

trang 

6 

c 

IN  THE  form 

trang 

7 

c 

TmETA=5G(YY.2)  KHERE  THETAiPHl/PHlO 

trang 

a 

c 

The  USER  Can  SPECIFY  THE  FUNCTION  SGtYYiZ) 

trang 

9 

c 

thf  user  must  define  as  Functions  of  (YYtn  the  following 

trang 

10 

c 

SG.  SGY.  SGZ.  SGYY,  SGY2 

trang 

11 

c 

SEE  USERS  MANUAL  FOR  RESTRICTIONS  ANO  INSTRUCTIONS 

trang 

12 

c 

M IS  THE  Index  for  The  tangential  plane 

trang 

13 

c 

J=1.2.3  IS  A LINE  INDEX  FOR  THANG  QUANTITIES 

trang 

14 

c 

trang 

15 

common  nc.mc.x ,pinf .dinf.phio. ioyaw.pi tPAO 

NEWCOM 

1 

COMMON  yZ(3) .YPHI (3) fC (25) .CZ (25) .CPHI ( 25 ) . R ( 20 , 25 ) 

NEWCOM 

2 

common  0(20.25) .p (20.25) .U(20t25) .V(20t25) ,w (20.25) .ASO(20.25) 

NEWCOM 

3 

common  CU(4.20.25) .CUP(4.20.25) 

NEWCOM 

4 

c 

• 

••  END  OF  rlaNX  common  ••• 

C03CSS 

32 

common  /CTRanG/  NS30.SGD(25) .SGY0(25) .SGYY0(25) 

NEWCOM 

6 

1 .gymdy.gyymoy.gyipoy.gyyipuy 

CTRANG 

3 

2 .MCP 

NEWCOM 

7 

common  /CBOOYZ  Z.BZZ.BPHPHI .BZPMI .TANCO.OELZ 

CHODY 

2 

1 .PHI (25) .B(25) .aZ(25) .BPHI (25) .COSPHI (25) .SlNPHl (25) 

C800Y 

3 

COMMON  /BLN02/  theta. OY. TG4(3) »TG5(3) *TG6(3) 

BLK02 

2 

1 .X(20).XZ(20.2).XH(20.2).XPh1(20.2).Y(25) 

BLK02 

3 

2 .TF4(20.2).TFA(20.2),TF7(20,2) 

BLK02 

4 

c 

trang 

IB 

c 

This  routine  is  VERSION  1 OF  TRANG  CORRESPONDING  EITHER 

TRANG 

19 

c 

TO  NO  CLUSTERING.  I.E..  SG(YY.Z)=YY 

trang 

20 

c 

OH  THE  USER  HAS  READ  IN  THE  PHI"S 

trang 

21 

c 

Trang 

22 

IF  (NSGD  .EO.  0)  GO  TO  50 

trang 

23 

c 

trang 

24 

c 

IHE  USER  READ  IN  THE  PHI"S 

trang 

25 

c 

trang 

26 

IF  (M  ,nE.  MCP)  GO  TO  30 

trang 

27 

SGY=GYIP0Y  t SGYY=GYY1P0Y  i SGZ=SGYZ=0. 

TRANG 

26 

GO  TO  2 

trang 

29 

30 

IF  (M  ,nE.  0)  GO  TO  35 

Trang 

30 

SGY=GYMOY  i SGYY=GYYMOY  S SGZ=S3YZ=0. 

trang 

31 

GO  TO  2 

Trang 

32 

35 

SG=SGO(M)  S SGY=SGYO(M|  i SGYY=SGYY0 (M)  i SGZ=SGYZ=0. 

trang 

33 

GO  TO  1 

trang 

34 

c 

trang 

35 

c 

CORRESPONDS  TO  NO  CLUSTERING 

trang 

36 

c 

TRANG 

37 

50 

Su=YY  i SGY=I.  S SGZ=SGYY=SGYZ=0. 

trang 

36 

c 

trang 

39 

c 

• 

••  The  following  statements  Should  appear  in  all  VERSIONS  ••• 

TRANG 

40 

1 

THtTA=SG 

trang 

41 

2 

YPHI (J( =1.0/ (PmIO®SGY) 

TRANG 

42 

3 

YZ ( J) =-S6Z/SGY 

trang 

43 

4 

Tu4 ( J) =SGY 

TRANG 

44 

5 

TG5(J)=SGYY/SGY 

trang 

45 

6 

TGb(J)=5GYZ/SGY 

TRANG 

46 

return 

trang 

47 

ENTRY  TRANGw 

TRANG 

4B 

IVERSON=l 

TRANG 

49 

WRITE  (6.3000)  IVERSON 

TRANG 

50 

IF  (NSGO  .EQ.  0)  WHITE  (6.3010) 

trang 

51 

IF  (NSGD  .NE.  0)  WRITE  (6.3020) 

TRANG 

52 

3000 

FORMAT ( 1H0.20X.WPROGRAM  TRANG* . 6X .• VERS  ION* . 1 4 ) 

TRANG 

53 

3010 

F0HMAT(11X,»E0UAL  SPACING  IN  TANGENTIAL  DIRECTION*) 

trang 

54 

3020 

F0HMAT(11X.*THE  PHI"S  WERE  READ  IN  BY  THE  USER*) 

TRANG 

55 

RETURN 

trang 

56 

END 

TRANG 

S7 

■R. 

7. 

Listing  of  Delivery  Version  of  TRANC. 
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10.2  Input  of  Uneven  Mesh 

In  this  approach,  the  user  can  select  an  uneven  mesh  spacing 
In  the  4),  or  r,  or  both  directions  by  inputting  to  the  program  the 
mesh  point  coordinates.  4 and/or  x where  x is  the  normalized  radial 
coordinate 


(h  and  c are  the  local  radius  of  the  body  and  bow  shock,  respectively). 
Note  that  in  this  method,  the  clustering  in  each  direction  is  independent 
of  the  other  variable  and  z. 

ilhen  an  uneven  spacing  is  desired  in  the  radial  direction,  NSFD  is 
read  in  as  the  total  number  of  points  in  the  radial  mesh  (^0)  (see, 
sec.  7.2).  The  values  x(N)  for  1 1 N - NSFI)  are  input  from  cards 
(c.f.,  sec.  7.3).  Note  that  it  is  required  that 

x(N)  < x(N+l)  N=l,2, • • • ,NSFD-1 

and 

x(l)  = 0.0,  x(NSFD)  = 1.0. 


If  a uniform  mesh  in  the  radial  direction  is  desired,  then  NSFD  is 
read  in  as  zero. 

When  an  uneven  spacing  is  desired  in  the  4>  direction,  NSGD  is 
read  in  as  the  total  number  of  points  in  the  <(>  mesh  (?^0)  (see,  sec.  7.2). 
The  values  <1(1!)  (in  degrees)  for  1 £ M - NSGD  are  input  from  cards 
(c.f.,  sec.  7.5).  Note  that  it  is  required  that 
(fi(M)  < 4i(M+l),  M=1,2,-*',NSGD 
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i 


and 


<(>(1)  = 0.0 


(J^(NSGD) 
Further,  for 


I 180.0,  for  the  syminetrlc  case  (PHI0=it) 

I 360.0,  for  the  nonsymmetrlc  case  (PHIO=27r). 
the  nonsymmetrlc  case  it  Is  also  required  that 


(f(NSGD-l)  = 360.  - (f(2). 


Our  experience  has  shown  that  in  order  to  obtain  accurate  results 
using  this  method,  the  data  points  must,  in  addition  to  the  above 
restrictions,  have  a certain  degree  of  smoothness.  We  have  fo..nd  that 
good  results  are  obtained  when  the  numerical  second  differences  of  the 
data  are  "smooth".  The  selection  of  mesh  points  is  primarily  by  trial 
and  error.  We  have  found  that  a workable  approach  is  to  use  a reasonable 
"eyeball"  choice  to  begin  with  and  then  smooth  the  data  using  a "mesh- 
preprocessor"  (see  Appendix  A) . Generally  the  preprocessor  must  be 
used  several  times  before  the  desired  smoothness  is  obtained.  This 
procedure  alters  the  initial  mesh  array  but  even  after  several  smoothing 
operations  the  final  mesh  will  retain  the  concentrations  in  regions 
where  originally  desired. 


11.  BENT  CONE  CALCUlJ^TIONS 

11.1  Applicable  Body  Shapes;  The  computational  approach  used  in 
the  program  is  applicable  to  a family  of  body  shapes  depicted  in  Fig.  8. 
The  body  shape  is  required  to  satisfy  the  following  restrictions: 

(1.)  The  complete  body  shape  must  be  continuous;  however, 

discontinuities  in  slopes  and/or  curvature  are  allowed. 
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(il.)  In  the  bent  cone  section,  the  body  must  be  a spherically 

blunted  right  circular  cone  with  respect  to  the  bend  axis. 

(Hi.)  The  body  shape  in  the  transition  and  aft-body  sections  can  be 
arbitrary  but  must  be  given  with  respect  to  the  aft-body  axis 
with  origin  as  indicated  in  Fig.  8. 

(iv.)  2„  < z 1 z (see  Fig.  8) 

2 en  n ! 

The  parameters  6,  = THETABN,  a = ALNS,  z = 2EN,  r = RN,  z = ZBAR, 
n n en  n 

= HN  (c.f.,  Fig.  8)  are  used  internally  in  the  code.  They  must  be 
defined  (in  ENTRY  BODYR,  see  sec.  12.1  of  Ref.  1)  either  by  direct 
input  or  in  terms  of  other  input  parameters.  Note  that  the  bent  nose 
section  geometry  given  in  the  delivery  version  of  BODY  is  a special  case  of 
the  above. 

11.2  Computational  Method  - The  general  approach  used  for  bent  nose 
calculations  is  to  perform  the  computation  along  various  Z-axes  (see 
Fig.  9).  Each  Z axis  is  parallel  to  the  aft-body  axis;  the  first  is 

the  axis  through  the  center  of  the  spherical  tip.  In  the  course  of  the 
calculation  in  the  bent  cone  region,  the  axis  is  shifted  down  automatically 
when  required.  The  rezonlng  required  when  the  axis  is  shifted  is  done 
in  subroutines  SHFAX  and  SAFAXD  (see  sec.  11.7  of  Ref.  1).  The  user 
can  control  the  shifting  of  the  axes  by  certain  input  parameters  (see 
the  next  paragraph  for  details).  Note,  that  the  axis  is  always  shifted  to 
the  aft-body  axis  before  the  calculation  enters  the  transition  region. 

11.3  Use : When  a bent  cone  calculation  is  to  be  performed  the 

input  variable  IBN  (see  sec.  7.2)  must  be  read  in  as  one  (1).  The  user 
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must  supply  the  body  geometry  consistent  with  section  11.1,  above.  In 
BODY  and  BODYR.  The  calculation  proceeds  automatically;  however,  the 
user  can  control  the  7.  axis  shifting  using  the  Input  variable  CENUF 
(see  sec.  1.2).  The  Z axis  Is  shifted  when  either  of  the  following 
criteria  are  met: 

(1.)  the  z axis  becomes  closer  than  CENUF  to  the  body  (see  Fig.  9) 

(11.)  the  axial  location  Is  within  0.5  of  ZEN  (see  Fig.  9). 

In  case  (1),  the  axis  Is  shifted  down  by  the  amount  ZMAXS  = CENUF*TAN(THETABN) 
but  never  past  the  aft-body  axis.  In  case  (1),  the  axis  Is  shifted  to 
the  aft-body  axis. 

Note  1:  I-Zhen  starting  the  calculation  from  spherical  blunt  body  data, 

care  should  be  exercised  so  that  the  starting  value  of  Z (l.e.,  the 
Initial  plane)  should  not  exceed  l.-SIN(THETABN+ALNS) . 

Note  2:  The  calculation  can  be  restarted  In  the  bent  cone  region  but 

= DELIT  (see  Fig.  9)  must  be  Inputted  (see  sec.  7.2). 
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12.  PROGRAM  STOP  .3 


The  normal  termination  of  the  program  is  when  either  z > ZEND  or 
K ^ KA.  The  program  can  be  prematurely  terminated  either  by  programmed 
stops  or  by  error  mode. 

12.1  Programmed  Stops:  The  various  programmed  stops  in  the  code 

are  listed  below  according  to  the  on-line  printed  message. 

(1)  NO  DATA  ON  TAPF.3 STOP 

(2)  THE  LAST  K ON  TAPE15  IS  (value  of  K)  LESS  THAN  (value  of  KSTART) 
Stops  1 and  2 indicate  that  the  vjrong  input  tape  is  being  used. 

(?)  DZ  IS  LESS  THAN  l.E-4 STOP— 

This  stop  is  caused  by  the  step  size  becoming  to  small.  Generally, 
it  indicates  that  the  axial  flow  component  is  approaching  the  sonic 
value. 

(4)  IN  SUBROUTINE  SHOCK  Cl  = (value  of  Cl)  L,M,K  = (values  of  L,M,K) 

This  message  is  printed  in  subroutine  SHOCK  when  a change  of  sign 
in  Cl  (see  sec.  3.3,  ref.  1)  has  been  detected.  Generally,  it  indicates 
either  subsonic  axial  flow  behind  the  boxj  shock  or  an  unrealistic  bow 
shock  (c.f.  Appendix  A,  ref.  1).  Here  L = 0 or  1 depending  on  whether 
the  calculation  is  in  the  predictor  or  corrector  step,  respectively. 

(5)  ITERATION  LIMIT  EXCEEDED 

(6)  LIMIT  ON  SHOCK  ITERATIONS  EXCEEDED 

Stops  5 and  6 are  executed  from  subroutine  DECODE  when  the  real  gas 
decoding  procedure  requires  more  than  LCNT  iterations  (5  is  printed  when 
this  occurs  at  an  interior  point;  (S  when  this  occurs  at  a boxj  shock 
point)  . 

(7)  OBLIQUE  SHOCK  ITERATION  EXCEEDS  LIMIT 

This  message  is  printed  in  subroutine  JLTfP  when  the  real  gas 
Iterations  for  the  obilque  shock  exceeds  2*LCNT. 

(8)  "EXIT  CALLED  ON  TENTH  FAILURE" 

This  message  is  printed  in  subroutine  RCAS  after  the  thermodynamic 
properties  are  out  of  range  10  times  (real  gas  calculations  only). 
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(9)  OUTSIDE  GAS  TABLES 

This  message  is  printed  in  subroutine  HRGAS  when  the  thermodynamic 
properties  are  out  of  range  (real  gas  only). 

Note:  When  the  messages  3 or  4 are  printed,  the  program  executes 

subroutine  SAVE  (see  sec.  12.4,  ref.  1)  before  the  run  is  terminated 
(this  does  not  require  inclusion  of  the  system  routine  RECOVR) . 
Therefore,  the  flow  field  is  printed  for  the  last  |IERRPR|  steps  before 
the  stop.  Also,  the  wall  pressure  and  aerodynamic  data  are  printed. 

12.2  Error  Mode:  At  various  locations  in  the  code,  errors  can 

cause  the  program  to  terminate.  These  errors  can  be  located  using  the 
system's  error  message  and  dump.  Many  times  the  cause  of  the  error 
can  be  determined  using  the  printout  obtained  from  the  recovery  routine 
RECOVR. 

A program  update  is  available  which  provides,  for  the  more  common 
occurring  errors,  programmed  stops  consisting  of  an  error  message  and 
a call  to  SAVE  (without  RECOVR,  see  the  note  in  the  previous  subsection) 
before  terminating  the  program.  This  update  is  given  in  Fig.  10. 

Most  of  the  stops  created  by  the  update  are  caused  by  a negative 
argument  in  SQRT.  In  these  cases,  the  error  messages  are  of  the  same 
general  form;  l.e.,  IN  ROUTINE  (subr.  name)  AT  CALL  NO.  (test  no.  within 
the  subr.)  NEGATIVE  SQRT  ROOT  of  (value  of  argument)  FOR  Z,K,M,N  (value 
of  z)  (station  no.)  (Y  mesh  index)  (X  mesh  index).  The  stops  with 
this  message  are: 
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•lOENT  DUMP 
•INSEMT  03CSb.lh0 

IERMPR=:IAMS(IER«PW) 

•BEFORE  EVAl.E-* 

IF  (OUHV  .GE.  0.)  GO  To  1001 

CALL  UPPbOPI ( AHE VAL • 1 • ? iR  tM.NtDUMl ) 

1001  CONIINOE 
•BEFORE  OECOOE.El 

IF  (CVU.l.M)  .GI.  -600.  .ANO.  CV(l.ltM)  ,LT,  700.)  00  TO  2001 
WRITE  (6.3*66)  M.CV(l.l.M) 

3*66  FORMAT ( IMl .•IN  SUBROUTINE  DECODE  THE  LOG  OF  PRESSURE  ON  PLANE^.I3. 

1 • ON  the  boot  IS^.1PE15.6.5X.^—  STOP  •) 

CALL  SAVE lOUM.UUM.DUM) 

2001  CONTINUE 
•BEFORE  OECOOE.SR 

IF  (NTFST  .LT.  0)  CALL  SAVE (UUM.OUM.DUM) 

•DELETE  DECODE. 26 

COUMP=QSQ»T3-U3«U3 

IF  (COUMP  .GE.  0.)  GO  TO  1001 

call  DmoSORT (6HUEC00E. 1 .2.K.H. 1 .COUMP) 

1001  WM=SQRT (COUMP) /a2 (J) 

•DELETE  0EC00E.7R 

COUMP=l.-CRAP2 

IF  (COUMP  .GE.  0.)  GO  TO  1002 

CALL  OMPSQRT (6H0EC00E.2.2.K.M.N.CDUMP) 

1002  CRAP3=CRAP2/( ( 1 ..SORT (COUMP) ) agEF) 

•delete  DECODE. 1 73. DECODE . 17* 

*8  CDUMP=CVN2«  (CVN2»0IT(E  2.P '•‘■'•DINF^  ( 2.  .8  , •GE^GO/GINl ) ) 

1 »PINF*PINF«GAMMA«GAMMA 

IF  (COUMP  .GE.  0.)  GO  TO  1003 
call  DMPSORT (6HOECOOE.3.Z.K.M.NC.COUMP) 

1003  PNM= (PINF .01NF»CVN2»S0Rr (COUMP) ) /(2.^GE ) 

•delete  decode. 1R2 

C0UMP=CVN2/CVTN 

IF  (COUMP  .g£.  0.)  GO  TO  100* 

CALL  DMPSORT (6H0EC00E t « . Z .K.N.NC.CDUMP ) 

100*  DN“=OINF»SORT (COUMP j 
•DELETE  DECODE. 216 

CDUMP=CVN2/CVT2 

IF  (COUMP  .GE.  0.)  GO  TO  1005 

CALL  DMPSORT (6MOECOOE.S.2.R.M.NC.COUMP) 

1005  DNM=OlNF»SORT (COUMP) 

•DELETE  RALL.AO 

COUMPi (E I A» ( 1 .♦BPH0B»»2) » («W*B2M) ••Z) /ASOH 

IF  (COUMP  .GE.  0.)  GO  TO  1001 

CALL  DMPSORT (*MRA|  L.l.Z.K.M.l.CDUMP) 

1001  BETA=SORT (COUMP) 

•delete  SMOCK. 20 

UV»US-CPM1C»VS 

COUMP = U Ta«cmu1.UV»*2) /ASQS 

IF  (COUMP  .OF.  0.)  GO  TO  loot 

CALC  DMPSORT (5HSM0CK. I . Z «R .M.NC .COUMP) 

1001  BET T*SURT (COUMPI 
•ADOFlLE . 

•DECK  DMPSORT 

SL'HROUTINE  OMPSORI  (NAME. RNTiZ.R.M.N. VALUE) 

C 

c name  is  The  name  of  The  routine 

C RNT  IS  IMF  NUMBER  FROM  WHICH  NAME  WAS  CALLED 

C Z IS  THE  Z value 

C R IS  THE  STATION  NO. 

C M IS  THE  PLANE  NO. 

C N IS  THE  RADIAL  POINT  NO. 

C VALUE  IS  THE  ARGUMENT  OF  SORT  ROOT 

C 

WRITE  (6.3000)  NAME.RNT.VAIUE.Z.R.H.N 
3000  FORMAT ( IHl .•IN  ROUTINE  •.AIO.^  AT  CALC  NO.^iI3,2X. 

1 •NEGATIVE  SORT  ROOT  OF^.IPE)5.6. 

2 • FOR  Z.R.M.N^. IPt 15.6.315) 

CALL  SAVE (DUM.OUM.OUM) 

STOP 

END 


0.  I.istinp  of  th(’  Frror-Modo  FpdaLe 
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Call 


Routine 

No. 

Description 

EVAL 

1 

in  evaluation  of  stability  condition 
(c.f.  sec.  3.6,  ref.  1),  indicates 
subsonic  axial  flow  or  neg.  sound  speed. 

DECODE 

1 

in  eq.  for  axial  velocity  component  at 
body  surface  (see,  sec.  3.4,  ref.  1), 
indicates  unrealistic  surface  flow 

DECODE 

2 

in  decoding  formula  for  interior  pts. 

(see  sec.  3.2,  ref.  1),  indicates  subsonic 
axial  flow 

DECODE 

3 

in  eq.  for  pressure  at  bow  shock  wave  (see 
sec.  3.3,  ref.  1),  can  only  occur  in  a real 
gas  run 

DECODE 

in  real  gas  iteration  formulas  at  bow  shock 
wave 

WALL 

1 

in  eq.  for  6i  (see  sec.  3.4,  ref.  1),  indicates 
subsonic  axial  flow  at  body  surface 

SHOCK 

1 

in  eq.  for  Bq  (see  sec.  3.3,  ref.  1),  indicates 
subsonic  axial  flow  at  bow  shock  wave 

Two  additional  stops  are  created  by  the  update.  For  one  the  following 
message  is  printed: 

IN  SUBROUTINE  DECODE  THE  LOG  OF  PRESSURE  ON  PLANE  (Y  index)  ON  THE 

BODY  IS  (value  of  log  of  pressure) STOP 

This  message  is  printed  in  DECODE  when  the  value  of  the  log  of  the  surface 
pressure  is  out  of  the  range  of  EXP.  The  other  stop,  also  executed  in 
DECODE,  can  occur  only  in  a real  gas  calculation.  The  printed  message  is 
P STAYS  NEGATIVE  AFTER  AVERAGING  Z,N,M,K  (value  of  z)  (X  mesh  index) 

(Y  mesh  index)  (station  no.) 

Note  that  this  message  can  appear  in  a perfect  gas  calculation  but  the 
program  is  not  terminated. 
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12.3  Remark: 

When  the  program  terminates  prematurely  it  generally  means  that 
either  an  error  has  been  committed  in  setting  up  the  run  or  a nonuniformity 
has  developed  in  the  flow  field  which  the  present  method  can  not  properly 
resolve.  In  either  case,  inaccurate  or  non-physical  results  are  present 
which  ultimately  produce  the  error  mode  or  programmed  stop.  Keep  in  mind, 
however,  that  it  is  possible  that  an  accurate  and  physically  consistent 
calculation  can  still  be  prematurely  terminated.  Obviously,  there  are  body 
shapes  and  free  stream  conditions  for  which  the  flow  field  will  violate 
the  basic  assumptions  of  the  code;  viz.,  that  the  flow  is  supersonic  in 
the  axial  direction  and  moreover,  that  it  remains  attached  to  the  body 
surface.  On  the  other  hand,  it  should  by  now  be  clear  that  not  every 
calculation  can  be  considered  accurate  and  physically  realistic  simply 
because  no  premature  termination  occurs.  All  results,  particularly  those 
for  new  classes  of  body  shapes,  should  be  carefully  studied  for  physical 
consistency.  Good  lucki 


U 

■1 


i 
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APPENDIX  A 

Mesh  Smoothing  Preprocessor 


Identification:  DIFF 

Important  Note:  DIFF  calls  two  subroutines  TRANGD  and  TRANFD.  These 

routines  bear  the  same  name  as  two  subroutines  in  D3CSS  because  they 
accomplish  the  same  thing.  But  note,  the  subroutine  structures  are 
different  and  thus  the  subroutines  should  be  kept  distinct. 
Description:  DIFF  Is  a separate  program  which  can  be  used  to  smooth 

the  mesh  coordinate  data  when  the  mesh  is  to  be  input  from  cards  (see 
sec.  10.2).  The  program  is  used  before  the  flow  field  calculation  and 
can  be  used  to  smooth  either  <j)  mesh  coordinates  or  x mesh  coordinates. 
The  basic  smoothing  operation  used  in  DIFF  is 


«l>..oothed  ■ <'*I«  * “l  * 

where  is  the  array  of  mesh  point  values  (either  i|>  (M)  or  x(N)) 

supplied.  The  smoothing  operation  retains  the  required  properties  of 
the  mesh  arrays  given  in  sec.  10.2.  When  x values  are  to  be  smoothed, 
the  last  value  of  x does  not  have  to  be  unity  since  DIFF  normalizes  the 
X values  after  input.  The  program  performs  any  number  of  such  smoothing 
operations  selected  by  the  user. 

Output : After  each  smoothing  C;peratlon,  DIFF  prints  on-line  the  mesh 

points  (or  x^)  and  the  numerical  approximations  to  the  first  and 
second  derivatives  of  the  underlying  mesh  clustering  functions  g^  and 
g^Y  fy  ^'’d  as  functions  of  the  mesh  Index  M (or  N) . Also, 
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after  the  last  smoothing  operation,  data  cards  are  punched  for  use 
in  the  flow  field  program. 

Input ; Input  to  DIFF  is  via  data  cards. 


Card  No. 

Format 

Variable 

Descrl 

ption 

1 

(15) 

N121 

No.  of  smoothing 

■ 

operations 

desired 

NXPHI 

No.  of  PHI 

(or  X 

to  be  read 

2 

(215) 

IXPHI 

=0,  PHI  are  read 

=1,  x's  are 

read 

3 

(5F10.10) 

XPHI (N) , 

the  PHI  (or  x ) 

N=l, NXPHI 

Sample  input  data  cards  for  the  three  typical  applications  of  DIFF 
are  given  on  page  45.  These  are:  smoothing  of  4>  for  the 
symmetric  problem,  smoothing  of  x , and  smoothing  of  4>  for  the 

nonsyrametrlc  problem,  respectively. 

Listings:  The  listings  for  DIFF  are  given  on  the  following  pages. 
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LISTING  OF  »*«E-PROCtSSON  PR0GM4K 


PROGHAM  DIFF ( INPUT  t OUT  PUT (PUNCH. TAPES* INPUT •TAPEbaOUTPUT) 
THIS  PROGRAM  DOES  l-^-l  Smoothing  of  data 

COMMON  /CTRAN/  NXPHI .XPHI 1200) .SGYO(200) .SGYYO(200) 

1 .GYMUY.GYYMDY.GYIPDY.GYYIPOY 

2 .MCP.PHIO.RAD.PI 

N121  IS  THE  NUMBER  OF  1-2-1  SMOOTHING  OF  DATA 
NXPHI  IS  THE  number  of  PHI"S  OR  X"S  TO  BE  READ 
ixPHi  = 0 Then  read  in  phi»s  (decrees) 

* 1 then  head  In  X"S 

note  that  the  X"S  do  not  have  to  be  read  in 

BETWEEN  0.  AND  I.  AS  THEY  ARE  NORMALIZED 

PI=*.»ATAN(1.)  » RAO=P1/180. 

S READ  (S.200U)  N121 
IF  )EOF(S))  999.10 
10  READ  (5.2000)  NXPHI. IXPHI 

READ  (5.ZI00)  (XPHI (N) .N=l .NXPhI ) 

2000  FORMAT) 1015) 

2100  FORMAT(5F10.0) 

PHlOO=XPHI (NXPHI ) S PHI0=PHI0D»RAD 
DO  25  N=1 .NXPHI 

25  XPHI (N) =XPhI (N) /XPHI (NXPHI ) 

N12IPI=N12I.l  S NKPHlKl~NXPHl-l 

DO  100  I=1.N121P1 

IM1=I-1 

WRITE  (b.ToOO)  IMl 

3000  formatuhi.'The  derivatives  after*.  13. 2X, 

1 *1-2-1  SMOOIHINGS  ARE  AS  follows*) 

IF  (IXPHI  ,E0.  0)  call  TRANGD 
IF  (IXPHi  .NE.  0)  CALL  TRANFD 
IF  (I  .EM.  N121P1)  GU  TO  125 
XPhINM1=xPHI ( 1 ) 

DO  50  N=2.NXPhIM1 
XPHIN=XPMI (N) 

XPHI (N) = (XPHINMI .2.*XPHIN»XPHI (N*l ) ) /♦, 

XPHINM 1 =XPH  IN 

50  continue 

IF  (IXPHI  .NE.  0)  GO  TO  100 
IF  (PHIO  .LE.  2.*PI-1.E-6)  go  TO  100 
XPHI (2) sXPHi (2) .360.-XPHI (NXPHIMl ) 

XPHI (NXPHIMI )=360.-XPHI (2) 

100  CONTINUE 

125  IF  (IXPHI  .NE.  0)  60  TO  200 
DO  150  N=l. NXPHI 
150  XPHI (N) *PhIOD*XPhI (N) 

200  PUNCH  AlOO.  (XPHI (N) .N=l .NXPHI I 
♦100  FORMAT (5F 10.5) 

GO  TO  5 
999  STOP 
END 
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J.ISIJMG  OF  ►'Hfc -PHOCESSOH  PPOCR/.M 


SUHHOUTINF.  TRANFO 
C 

c TPANFD  defines  quantities  needed  bv  subroutine 

C TRANF  KHtN  THE  USER  READS  IN  THE  SFIX.YtZ)  DATA  POINTS 

c 

COMMON  /CTRAN  / NSFO.SFD  (200)  tSr>iD(200l  tSFAXOiaOOl 

c 

NSFDM1=NSF0-1 
DX=1 ./Float (NSFDMl ) 

TNOOX=:l  ./ (2.»DX)  t DXSQ=I  ./DX*»2 
DO  SO  N=2. NSFDMl 

SFXO(N)=(SrO(N«l» -SFU (N~I ) ) ‘TwOOX 
SO  SFXXD  (N)  = tSFO  (N*  1 ) -2.  «SFO  (N)  ♦SFD  (N-U  ) ‘OXSO 

c 

C SFXXD  IS  ASSUMED  LINEAR  ON  (0.20X)  AND  Cl-2DX,1) 

C 

SFXXDU)=2.*SFXXO(2)-SFXXD(3) 

SFXXO(NSFD) =2. »SFXXD (NSFDMl ) -SFXXD (NSFD-2) 

SFXD( 1 ) *SFXD(2)-.S»DX» (SFXXDI I) ♦SFXXD(2) ) 

SFX0(NSFD)=SFX0(NSFDM1) ♦.5»DX«(SFXX0(NSFD)*SFXXDINSEDM1I) 
WRITE  (6>J300) 

3300  FORMAT ( lHa.AX,«N»t 15X,»SF». I7X,»SFX»t 16Xi»SFXX») 

DO  125  N=: ,NSFO 

125  WRITE  (0.3X00)  N . SFO (N ), SF XO (N) .SFXXD (N) 

3X00  FORMATdH  .I5.IPXE20.6) 

RETURN 

ENO 
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LISTING  OF  PWF -PHOCESSOM  PNOGHAM  (CONTINUED) 


SUBWOUTlNt  TBANGO 

tmango  defines  quantities  neeoeo  by  subroutine 
trang  »hen  the  user  heads  in  the  phi  values 

COMMON  /CTHAN  / NSGO.SGD (200) tSGYO(?00) tSGYY0(200) 

1 tGYMOY .GYYMOT.GYIPDYfGYYlPOr 

2 tMCP.PMlO.RAO.Pl 

NSc.  .l«NSGO-l 
Or=l ./FLOAT (NSGDMl ) 

TV00Y»1./(2.«0Y)  S 0YSQ=1./0V»«2 
DO  SO  M»2, NSGDMl 

SGYO (M) » (SGUiM. 1 ) -SGUIM-1 ) ) •TVOOY 

so  SGYY0(M)= (S&D(M*1)-2.»SGU(M)»SG0(M-I) )»DYSO 
IF  (PHIO  .Lt.  2.»(M-l.E-6)  GO  TO  7S 

note  that  SGD)Y«1)=SG0(T)»1  FOR  NON-SYMMETRIC  PROBLEM  (PHlO-360) 

SGYD( 1 ) =SGYO(NSGO) =(SGO(2)-SGO(NSGO-I ) ♦I.)»TV00Y 

SGYYDU  ) =SGYYU  (NSGD)  = (SGO(2)-2.»SGO(l)  •SGO(NS60-I)-I.)»OYSO 

GO  TO  100 

note  that  for  symmethic  problem  (PHio^ieo) 

SGYYD  IS  ASSUMED  linear  ON  (-0Y.2DY)  AND  (1-20Y*1«0Y) 

TS  SGYYOil)=2.»SGYYO(?)-SGYYO(3) 

SGYYO(NSGO)=2.*SGYYD(NSGOM1 ) -SGYYD (NSGO-2) 

GYYMOY=2. 'SOYTOI 1 )-SGYY0(2) 

GYY lPDY=^.•SGYYO (N5GD) -SGYYD(NSG0M1 ) 

SGYOU  ) =SGY0(2)--S»DY»  (SGYYOd)  ♦SGYY0(2)  ) 

SGYO (NSGO) =SGYO( NSGDMl ) •. 5»0Y» (SGYYD (NSGO) *S6Y YD (NSOOMl ) ) 
OYMOY=SGYO (2)-2.*DY»SGYY0( I ) 

GY1POY=SGYO (NSGOMl ) ♦ 2 . »D Y»SGYYO ( N5G0 ) 

100  WRITE  (S.IJOO) 

3300  format ( IMO. .T.»M»,15A.»PHI»»18X.»SG»,17X»»S6Y». JbXf •S6YY») 

M=0  $ WRITE  (6,3350)  M,GYMDYe('YYMDY 
33S0  FORMATdH  , IS.AOX,  1P3E20.6) 

PMrOO=PHIO/HAO 
DO  I2S  M*l,NSGO 
PMIM=SG0(M)»PHI00 

12S  WRITE  (6,3400)  M ,PH IM , $60 ( M) ,SGYO (M) .SGYYD (H) 

3400  FORMATdH  .15,1P4E20,6) 

M=NSGO*l  i WRITE  (6,3350)  N.GYIPOY.GYYIPOY 

return 

END 


I 
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LISTING  OF  QAT*  CAROS  FOR  PRE-PHOCESSOR  PROORAN 


10 

NUHUER 

OF 

1-2-1 

SMOOTHINGS 

33 

0 

NUMBER 

OF 

PMI"S 

TO  READ  ; 

0. 

A. 

8. 

12. 

16.  j 

20. 

2A. 

28. 

32. 

36.  1 

AO. 

AA. 

A8. 

52. 

56.  J 

AO. 

6A. 

68. 

72. 

76. 

no. 

8A. 

88. 

92. 

96.  j 

100. 

lOA. 

110. 

119. 

132.  1 

IA8. 

10 

16A. 

180. 

NUMBER 

OF 

1-2-1 

SMOOTHINGS  j 

26 

1 

NUMBER 

OF 

X**5  TO  READ  i 

0. 

1. 

2. 

3. 

A. 

5. 

6. 

7. 

8. 

9. 

10. 

12. 

lA. 

16. 

19. 

22. 

26. 

30. 

35. 

AO. 

AS. 

SO. 

S5. 

60. 

65. 

70. 

10 

NUMBER 

OF 

1-2-1 

SMOOTHINGS 

65 

0 

number 

OF 

PHI»S 

TO  READ 

0. 

A. 

8. 

12. 

16. 

20. 

2A. 

28. 

32. 

36. 

AO. 

AA. 

A8. 

52. 

56. 

60. 

6A. 

68  . 

72. 

76. 

no. 

8a. 

88. 

92. 

96. 

100. 

lOA. 

110. 

119. 

132. 

1A8. 

16A. 

180. 

196. 

212. 

228. 

2A1. 

2S0. 

256. 

260. 

26A. 

268. 

272, 

276. 

280. 

28A. 

288. 

292. 

296. 

300. 

30A. 

308. 

312. 

316. 

320. 

32A. 

328. 

332, 

336. 

3A0. 

3AA. 

3A8. 

352. 

356. 

360. 
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APPENDIX  B 
Sample  Run 

In  this  Appendix,  a sample  calculation  is  presented  to  illustrate 

the  use  of  the  code  for  a realistic  reentry  configuration.  The  body 

shape  considered  is  a spherically  blunted  20“  - 5“  biconic  with  a 

windside  flat  and  flap  (see  Fig.  B-1) . A perfect  gas  is  assumed  with 

y = 1.4.  The  freestream  conditions  are: 

= 15  (Mach  number) 

2 

p^  = 1.0  Ibj/ft  (pressure) 
p = l.xlO”^  (lb,-sec^) /f t^  (density) 

GO  f 

a = 7.5“  (angle  of  attack) 

6=0“  (side  slip  angle) 

The  blunt  body  region  was  determined  using  the  program  of  ref.  2.  The 

initial  plane  was  placed  at  the  sphere-cone  juncture  (z  = .658). 

In  setting  up  this  particular  calculation,  it  is  desirable  to  use 

the  following  computational  options: 

i.)  wall  entropy  reduction  (c.f.,  sec.  4.2,  ref.  1) 

ii.)  the  Mod  3 version  of  the  wall  boundary  conditions  with 
second  order  accuracy  (c.f.  sec.  3.4  and  10.7,  ref.  1) 

ill.)  a meridional  mesh  spacing  of  = 4“  in  the  vicinity  of  the 
flat  and  flap  but  a coarser  mesh  elsewhere 

The  other  parameters  and  options  used  in  the  calculation  can  be  found 

from  the  enclosed  input  and  output  which  follows.  Note  that  items  i.) 

and  il.)  are  subject  to  Internal  program  controls  when  the  body  has 

discontinuous  slopes  as  in  the  present  example  (c.f.  sec.  4.1,  ref.  1) . 


-46-  I 

J 


NSWC/WOL/TR  77-32 


Fig.  B-1,  Body  Geometry  for  Sample  Run 
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This  calculation  is  performed  in  three  separate  runs;  viz.. 

Run  No.  1 (sphere-cone  juncture  _ z 1 20'*'):  In  this  run,  the  meridional 

mesh  is  uniform  with  A(p  = 15°.  Note  that  after  the  biconic  juncture 
(z  = 15),  the  option  ii.)  is  automatically  modified  in  the  code. 

Run  No.  2 (20  f z - 33  ) : The  purpose  of  stopping  and  restarting  the 
calculation  near  z = 20  is  twofold.  First,  the  wall  point  calculation 
procedure  is  changed  back  to  the  Mod  3 version  with  second  order 
' accuracy.  It  was  felt  that  z = 20  is  sufficiently  far  downstream  of 

I 

j the  biconic  juncture  for  this  to  be  done  "safely".  The  other  purpose 

I 

is  to  perform  a preliminary  rezone  to  a uniform  meridional  mesh  with 
A4^  = 7.5°.  This  is  done  in  anticipation  of  the  final  reduction  of  the 
mesh  to  Atj>  = A°.  Generally  it  is  safer  computational  practice,  when 
rezoning  for  a finer  mesh,  to ^perform  this  in  stages  so  that  in  each 
rezone  the  mesh  lengths  in  each  direction  are  reduced  by  not  more  than 
one-half  their  previous  lengths. 

Run  No.  3 (33*'  f z f 60):  In  this  run,  the  mesh  is  rezoned  for  A<(>  = 4° 

near  the  flat  and  flap.  The  meridional  mesh  used  is  nonuniform  and 

I 

: input  from  data  cards  (c.f.,  sec.  10.2)..  It  has  A<j>  = 4°  for  0 £ <1)  - 92° 

I 

with  a gradual  increase  to  Aifi  = 16°  at  <}>  = 180°. 

The  input  data  cards  and  the  printed  output  (abbreviated)  for  each 
run  are  given  in  the  following  pages. 

Note  that  the  run  is  terminated  and  subsequently  restarted  at  the 
first  marching  step  beyond  this  value  of  z. 
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